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ABSTRACT 

We present the results of a detailed analysis of 121 ground-based high-resolution high 
S/N spectroscopic measurements spread over 3 years for the (3 Cephei star 9 Ophiuchi. 
We discovered 9 Oph to be a triple system. In addition to the already known Speckle 
B5 companion of the B2 primary, we showed the presence of a low-mass spectroscopic 
companion and we derived an orbital period of 56.71 days with an eccentricity of 
0.1670. After removing the orbit we determined two frequencies for the primary in the 
residual radial velocities: f\ = 7.1160 cd _1 and fi = 7.4676 cd~ . We also found the 
presence of fs — 7.3696 cd _1 by means of a two dimensional frequency search across 
the Si III 4567 A profiles. We identified the m- value of the main mode with frequency 
/i by taking into account the photometric identifications of the degrees £. By means 
of the moment method and the amplitude and phase variations across the line profile, 
we derived (4,mi) = (2,-1). This result allows us to fix the mode identifications 
of the whole quintuplet for which three components were detected in photometry. 
This is of particular use for our forthcoming seismic modelling of the primary. We also 
determined stellar parameters of the primary by non-local thermodynamic equilibrium 
hydrogen, helium and silicon line profile fitting and we obtained T g = 24000 K and 
log<? = 4.1, which is consistent with photometrically determined values. 

Key words: stars: variables: other - stars: early-type - stars: oscillations - stars: 
individual: 9 Oph - techniques: spectroscopic 



1 INTRODUCTION 

Recently (3 Cephei pulsators have proven to be very promis- 
ing targets for asteroseismic purposes. Besides the derivation 
of strong constraints on the mass and metallicity, among 
other global parameters, the fitting of three independent 
observed frequencies constrained the core overshooting pa- 
rameter of the j3 Cephei stars V836 Cen and v Eri. More- 
over, from the observation of two multiplets, non-rigid ro- 
tatio n was proven for both stars (Aerts et al. l2003l . Dupret 
et al. l2004l and Pamyatnykh et al.|2004) . In addition, it has 
been shown that non-standard stellar models have to be in- 
voked in order t o fit four independent frequencies of v Eri 
(Ausseloos et al. l2004l) . 

* Based on observations obtained with the CORALIE, FEROS 
and GIRAFFE echelle spectrographs attached respectively to the 
1.2m Leonard Euler telescope (La Silla, Chile), to the ESO 2.2m 
telescope (La Silla, Chile) and to the SAAO 1.9m telescope (South 
Africa) . 

f E-mail: maryline@ster.kuleuven.ac.be 
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By increasing the number of such asteroseismic studies 
of additional /3 Cephei target stars, it is clear that our un- 
derstanding of the internal structure of these massive stars 
will significantly improve in the near future. In this frame- 
work a three-site photometric campaign as well as a long- 
term spectroscopic monitoring were dedicated to the B 2 IV 
(3 Cephei star 9 Ophiuchi (HD 157056, V ma9 = 3.248). The 
detailed analyses of the gathered photometric and spectro- 
sco pic obs ervations are respectively described by Handler et 
al.| 1120051) . hereafter Paper I, and in the present paper. 

In Paper I, the frequency analysis and mode identifica- 
tion from 1303 colour photometric measurements show the 
presence of one radial mode, one rotationally split £ — 1 
triplet and possibly three components of a rotationally split 
I — 2 quintuplet. Such a pulsation spec trum is similar to the 
one observed for V836 Cen (Aerts et al. l2004l) . Consequently 
fruitful results can be expected from our asteroseismic study 
devoted to 6 Oph. The present paper describes the detailed 
analysis of our spectroscopic data with the main aim of iden- 
tifying the m- value of the mode with the highest amplitude, 
belonging to the quintuplet. In that way we can add an in- 
dependent well identified mode to the radial one and the 
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central mode of the triplet already detected. Such an iden- 
tification will allow seismic modelling of the star, which will 
be done in the near future. 

The spectroscopic variability of 9 Oph, as well as the 
main period of some 0.14 days, have been known for a long 
time but no spectroscopic mode identification was yet at- 
tempted. Radial velocity variations with a m ain p eriod of 
0.2862 d were first discovered by Henroteau (|l922j). From 
new radial velocity measurements, Van Hoof et al.| <ll956l) 
suggested a period of 0.15 d, which corresponds to about half 
Henroteau's period. They also noticed that the mean veloc- 
ity at the time of their observations was —8 km s" 1 , a value 
distinctly different from Henroteau's values of 1920 and 
1922, which were and —15 km s _1 , respectively. McNa- 
mara| (Il957f) and afterwards Van Hoof and Blaauw| |l95St) 
confirmed the period found by Van Hoof et al.| jl956T) and 
refined it to 0.1404 d. 

From their radial velocities measurem ents o f members 
of the Sco-Cen association Levato et al.| lll987l) suspected 
6 Oph to be a single-lined spectroscopic binary but they 
could not derive orbital elements. Moreover, McAlister et al.| 
( 1993) detec ted 9 Oph to be a Speckle binary and Handler 
et al.| J2005T) determined that the Speckle companion must 
be a B5 main sequence star. 

In the literature, values for the projected rotational ve- 
locity v sin i of 8 Oph were derived. Brown and Verschuerenl 
1 1997) measured vsini = 31 ± 3 km s _1 and Abt et al.| 
i 2002) found vsini ~ 30 km s _1 for our studied star, corn- 



Table 1. Observing logbook of our spectroscopy of Oph 



Instrument Number of JD 

observations 2450000 + 
Start End 



patible with the previous result. 

The paper is organized as follows. In Section 2 we give 
a description of our observations and data reductions. In 
Section 3 we show that the primary has a low-mass spec- 
troscopic companion besides the distant Speckle companion 
and we derive orbital parameters. Section 4 is devoted to 
the frequency analysis on our data after removing the orbit, 
followed by the identification of the wavenumbers (£, m) for 
the main frequency /i = 7.1160 cd _1 . In Section 5 we derive 
the position of 6 Oph in the Hertzsprung-Russel diagram. In 
Section 6 we end the paper with a summary. 



2 OBSERVATIONS AND DATA REDUCTIONS 

Our spectroscopic data were obtained with the CORALIE 
echelle spectrograph attached to the 1.2m Leonard Euler 
telescope in La Silla in Chile, during several runs spread over 
2000-2003. We also have at our disposal some measurements 
gathered with the FEROS and GIRAFFE echelle spectro- 
graphs attached respectively to the ESO 2.2m telescope (La 
Silla, Chile) and to the SAAO 1.9m telescope (South Africa). 
The number of observations and the ranges of their Julian 
Dates are given in TableEl The two last runs are respectively 
a FEROS run and a GIRAFFE run. 

An on-line reduction of the CORALIE spectra, using 
the INTER- TACOS software package, is available. For a de- 
scription of the reduction process we refer to Baranne et 
al.| lll99 6). We did a more precise correction for the pixel- 
to-pixel sensitivity variations by using all available flatfields 
obtained during the night instead of using only one flatfield, 
as is done by the on-line reduction procedure. The FEROS 
data reduction was performed using the on-line FEROS re- 
duction pipeline which makes use of the ESO-MIDAS soft- 
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ware package. This reduction includes consecutively the ex- 
traction of the inter-order background of the echelle spectra, 
the straightening and extracting of the orders, the removal 
of the blaze function and the pixel-to-pixel variations, the re- 
binning of the orders to the wavelength scale by using ThAr 
calibration spectra and a merging of the echelle orders. We 
performed an additional correction for the wavelength sen- 
sitivity of the shape of the internal flatfields by means of a 
smoothed average of dome flatfields. The GIRAFFE spec- 
tra were reduced by using the GIRAFFE pipeline reduction 
program XSPEC2, which follows a similar procedure as the 
FEROS pipeline reduction. The spectra of a ThAr arc lamp, 
taken at regular intervals to calibrate possible drifts in wave- 
length, and two different types of flatfields were used for 
the reduction. We note that bias frames were not separately 
taken as the bias value was calculated from the overscan 
region on the CCD. Flatfielding was accomplished by using 
camera flats, which were obtained by illuminating the cam- 
era with uniform light using a tungsten filament lamp and 
a diffusing screen. The blaze correction was determined by 
measuring the response across each order when the fibre was 
illuminated by the tungsten lamp (fibre flats). 

Finally, all spectra were normalized to the continuum 
by a cubic spline function, and the heliocentric corrections 
were computed. 



3 A TRIPLE SYSTEM 

For our study of the line-profile variability we used the Si III 
triplet around 4567 A. We point out that an effective tem- 
perature of some 18400 K was derived in Paper I for the 
B5 Speckle companion, which implies that this companion 
contributes very little to the Si III lines. Mor eover , the as- 
trometric data present ed by McAlister et al.| Jl993l) and by 
Shatsky & Tokovinin| i2002f) point towards an orbital pe- 
riod of order 100 yr, so our approximation that it shows 
no radial-velocity variation on a timescale of 3 years seems 
justified. 

We computed the first three velocity moments < v 1 >, 



< v 1 > and < v 3 > (see Aerts et al.| (11992ft for a definition 
of the moments of a line profile) of the Si III 4553 A line with 
the aim of performing a frequency analysis. The integration 
boundaries of the moments were dynamically determined by 
visual inspection of each spectrum in order to avoid noisy 
continuum. To perform the frequency analysis on the first 
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Figure 1. Phase diagram of the radial velocity placed at average zero with respect to the orbital period P or b = 56.71 d. The line 
corresponds to the orbital solution obtained with FOTEL (Tabled! . 



moment < v 1 >, which is the radial velocity placed at av - 
erage zero, we used the PDM metho d fStellingwer fll97c?) 
and the Scargle method (Scargle ll98ll) . We tested frequen- 
cies from to 10 cycles per day (c d _1 ) with a frequency 
step of 0.0001 c d _1 . The error estimate of our determined 
frequencies is between 0.0001 c d _1 and 0.001 c d _1 . We 
obtained the same results with both methods. 

The frequency that dominates the radial velocity is 
0.0175 c d _1 . Such a long period of some 57 days does 
not belong to the pulsation period range of /3 Cephei stars 
and is too long to be the rotation period of a star with 
vsini ~ 30 km s _1 and R ~ 5 Rq. Consequently we at- 
tributed this variability with a peak-to-peak amplitude of 
~ 25 km s _1 to binarity. 

In order to derive the orbital parameters of the binary 
system we used the code FOTEL developed by Hadrava| 
(1990). The program FOTEL is designed to solve either si- 
multaneously or separately light- and radial velocity curves 
of binary stars with a possible distant third component. For 
8 Oph we only used radial velocity measurements because 
there is no sign of the binarity in photometric measurements 
at our disposal in Paper I. We found an orbital period of 
Porb = 56.71 days with eccentricity e = 0.1670. The other 
orbital elements and their standard errors are listed in Ta- 
blc|21 The phase diagram with respect to P or b is shown in 

Figjn 

If we consider that the orbital inclination i or b corre- 
sponds to the inclination of the primary i rot , the result of 
the mode identification described in Section 4.2 leads to i or b 
in the interval [70°, 90°]. With a mass for the primary com- 
ponent of some 9 Mq (see Paper I), we conclude that the 
secondary component has a mass lower than 1 Mq . We note 
that the duty cycle of the photometric data in Paper I is not 
sufficient to detect a possible eclipse. 

We therefore conclude that 6 Oph is a triple system 
composed of a B2 primary, a low-mass spectroscopic sec- 
ondary and a Speckle B5 star. 



4 A PULSATING PRIMARY 

While the tertiary contributes significantly to the spectrum 
of the triple system, its Si III lines are so weak that they 



Table 2. The orbital parameters for 9 Oph. P or b is the orbital 
period, T the time of periastron passage, e the eccentricity, K\ 
the semi-amplitude of the radial velocity curve of component 1, 
u) the periastron longitude, v 7 the system velocity. 



Porb (days) 


56.712 ± 0.046 


To (HJD) 


51816.99 ± 2.71 


e 


0.1670 ± 0.043 


Ki (km s" 1 ) 


11.345 ± 1.099 


lo (degrees) 


128.046 ± 0.319 


Vj (km s" 1 ) 


7.202 ± 0.369 


ai sin i or b (a.u.) 


0.058311 


/i(m) (M Q ) 


0.0082243 


rms (km s" 1 ) 


3.35 



must be situated completely within those of the primary. 
We have no information about possible variability of the ter- 
tiary. In any case, should it be a rotational or a pulsational 
variable star, we expect its frequency range to be vastly dif- 
ferent from the one of the primary, given its mid-B spectra 
type. Its presence will therefore not affect the pulsational 
interpretation of the primary we present below. A similar- 
situation of having a companion spectral line within the one 
of th e prim ary occurs for the /3 Cephei star j3 Crucis (Aerts 
et a!. ll99St> . 



4.1 Frequency analysis 

After removing the orbit we searched for additional fre- 
quencies in the residual radial velocities by subsequent 
prewhitening. We found /i = 7.1160 cd _1 and fa = 
7.4676 cd _1 . Phase diagrams for the first moment are shown 
in Fig.|21 Note that /i is the frequency already detected 
in previous spectroscopic observations of 6 Oph and both 
frequencies correspond to those found in photometric mea- 
surements in Paper I. In the radial velocity and in the higher 
order moments we could not find any additional frequency 
after prewhitening with /i and fa. 

We then attempted to determine other frequencies by 
means of a two dimensional Scargle frequency analysis on 
the spectra themselves. After prewhitening with the main 
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Figure 2. Upper panel: phase diagram of the first moment com- 
puted from the Si III 4553 A line, after removing the orbit, 
for /i=7.1160 c d — 1 . Lower panel: phase diagram of the data 
prewhitened with /i for fa = 7.4676 c d — 1 . 



frequency fa clearly present in the data, we found the pres- 
ence of both frequencies fa and fa = 7.3696 cd -1 . The latter 
was also observed in photometry of Paper I. The amplitude 
and phase variations across the Si III 4553 A line are shown 
in Fig.EJfor fa, fa and fa. 



varied so that the theoretical amplitudes of the first moment 
correspond to the observed ones within the errors. 

4-2.1 The moment method 

We first used the moment method. With this method, the 
wavenumbers (£, m) and the other continuous velocity pa- 
rameters are determined in such a way that the theoretically 
computed first three moment variations of a line profil e best 
fit the observed ones. We refer to Briquet & Aerts| ll2003h 
for the latest version of the technique, which was improved 
by these authors in order to perform a simultaneous identi- 
fication of all the modes that are present in the data. Since 
we only detected fa = 7.1160 cd" 1 and fa = 7.4676 cd" 1 
in the first three moments and since fa is unambiguously 
identified as a radial mode by photometry (Paper I), we per- 
formed a multiperiodic mode identification with these two 
latter frequencies, the mode with frequency fa being fixed to 
(£2, Tri'i) = (0, ). According to the criterion given by Briquet 
& Aerts (2003), we only needed to test l\ ^ 3. 

The results of the mode identification are given in Ta- 
ble^ The lower the discriminant E, the better the agree- 
ment between theoretical and observed moment values. We 
can conclude that the mode identification by means of the 
moment method favours an m = —1. Since the photo- 
metric mode identification described in Paper I excludes 
an i\ — 3, the mode identification gives a preference to 
(4, mi) = (2, -1) or (1, -1). However, the solution (2, -2) 
remains possible. For the solutions in italic in Table^ we 
compare the moment values with the observed ones in Fig.^] 
in order to check if we can choose between these three pos- 
sibilities. We conclude that we cannot discriminate between 
(d,mi) = (2, —1) or (1, —1) and the solution (2, —2) gives a 
theoretical second moment amplitude slightly too large com- 
pared to the observed one. We also computed observed and 
theoretical moments up to order 6 in order to attempt to dis- 
criminate between the different solutions, as could b e don e 
for several slowly pulsating B stars by De Cat et al. (2005). 
The solutions with (2, —1) and (1, —1) give very similar first 
six moments and the solution with (2,-2) shows a slightly 
less good fit to the even moments. 



4.2 Mode identification 

In Paper I three frequencies of a rotationally split I — 2 
quintuplet were observed for 9 Oph, which are 7.1160 cd -1 , 
7.2881 cd" 1 and 7.3697 cd -1 . The former and latter corre- 
spond to fa and fa. Consequently the main I — 2 mode with 
frequency 7.1160 c d _1 has an m-value equal to —2 or —1. By 
means of our spectroscopic data we will lift this ambiguity. 

To perform our mode identification by different meth- 
ods we fixed certain parameters while others are fully free. 
We took a linear lim b-dark ening coefficient u of 0.32 (see 
e.g. Wade & Rucinski ll98fJ) . We used the if -value given by 
K = GM/uj 2 R 3 , where M is the mass, R the radius and 
u> the angular pulsation frequency. Note that slightly dif- 
ferent values of these parameters within the errors do not 
change the result of the mode identification. In order to de- 
fine a fine grid of unknown stellar and pulsation parameters 
without testing useless cases, the velocity amplitude A p was 



4-2.2 The amplitude and phase variations across the line 
profile 

In order to discriminate between the best moment solutions 
we also used the amplitude and phas e dist ributions shown 
in Fig.El (see Telting and Schrijvers| lll997l) . Schrijvers and 
Telting (1999) for a definition). Because of the zero ampli- 
tude at the center of the line and a corresponding phase shift 
of 7T, one can directly conclude that the mode correspond- 
ing to fa is axisymmetric. It corroborates the photometric 
mode identification. Moreover, because of opposite slopes for 
the phase distribution, the modes corresponding to fa and 
fa have m-values of opposite sign. We point out that this 
conclusion can be made for /3 Cephei stars having p-mode 
pulsations while it c annot be concluded for SPBs with g- 
modes (De Cat et al. l2005l) . 

In Fig.|3]we also compare the theoretical and observed 
amplitude and phase variations across the line profile for 
fa for the best solutions of the moment method given 
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Table 3. The ten best solutions of the mode identification through the discriminant £ (see Briquet & Aerts| 12003?) ) . using the Si III 
4553 A line. We took K\ = 0.120 and K2 = 0.109 for the ratios of the amplitude of the horizontal and of the vertical motion. A p is the 
amplitude of the radial part of the pulsation velocity, expressed in kms -1 ; fr.max and Dt.max are respectively the maximum radial and 
tangential surface velocity due to the two modes, expressed in kms -1 ; i ro t is the inclination angle, expressed in degrees; vsvai is the 
projected rotational velocity, expressed in kms -1 and a is the intrinsic line-profile width, also expressed in kms -1 . 
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Figure 3. Amplitude and phase distributions for f±, fi and for the Si III 4553 A line (dots) and comparisons for /1 with the theoretical 
ones for the solutions given in italic in TablelSI The solutions with (<i,mi) = (2, —1), (1, —1) and (2, —2) are respectively represented 
by a full line, a dashed line and a dot-dashed line. The amplitudes are unitless and the phases are expressed in it radians. 



m italic in Tabled] This comparison allows us to exclude 
(£i,mi) — (2,-2) since the theoretical variations do not 
mimic the observed ones for this solution. 

As an additional check we also performed an indepen- 
dent mode identification b y mea ns of the following proce- 
dure. By using Townsend's| jl997l) codes, called BRUCE and 
KYLIE, we generated theoretical line profile time series for 
several wavenumbers (I, m) and for a large grid of other con- 
tinuous parameters. Subsequently we computed the ampli- 
tude and phase variations across the line profile for each gen- 
erated line profile time series. Finally, by comparing these 
theoretically computed amplitude and phase distributions 
with the observed ones, we derived the couple (£, m) and 
other parameters which led to the best fit. 

Because generating line-profile variations is very time 



consuming, even with nowadays computers, we performed 
this procedure for the frequency /1 and for couples (2, — 1), 
(1, —1) and (2, —2) since our goal is to discriminate between 
these possible (I, m) for the main mode. We covered the 
parameter space by varying the free parameters in the fol- 
lowing way: the projected rotational velocity usini from 28 
to 32 kms -1 with a step 1 kms -1 , the inclination of the star 
i ro t from 1° to 90° with a step 1°, the line-profile width a 
from 1 to 20 kms -1 with a step 1 kms -1 . For each tested 
(£,m,i), we varied the velocity amplitude A p from 0.8 A p 
to 1.2 A* with a step 0.05 A p where A* is the value of the 
velocity amplitude which leads to a theoretical first moment 
amplitude equal to the observed one. 

In Fig.0 we show the best fit model compared to obser- 
vation. For a (2, —1) mode, one has vsini = 28 kms -1 , i rot 
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Figure 5. Comparisons between the observed amplitude distri- 
bution (dots) for /i computed for the Si III 4553 A line (dots) and 
the theoretically computed one for (li,m\) = (2, —1) (full line), 
for (h, mi) = (1,-1) (dashed line) and for (fi,rai) = (2,-2) 
(dot-dashed line). 



are (£i,mi) = (2,-1) and consequently the frequencies 
fi = 7.1160 ccT 1 , 7.2881 ccT 1 and f 3 = 7.3697 cd" 1 
of the rotationally split I = 2 quintuplet correspond to 
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Figure 4. Comparisons between the observed first three moments 
(dots) and three of the best solutions given in italic in Tablelsl 
The solutions with {t\,m\) = (2,-1), (1,-1) and (2,-2) are 
respectively represented by a full line, a dashed line and a dot- 
dashed line. 



= 78°, a =19 kms" 1 and A 1 = 31.62 kms" 1 . For a (1, -1) 
mode, the best solution is obtained for usini = 29 kms , 
irot = 70°, a = 19 kms -1 and A 1 = 10.20 kms -1 . Finally, 
for a (2,-2) mode, the parameter values are vsini = 32 
kms" 1 , i rot = 90°, a =19 kms" 1 and A 1 = 19.29 kms -1 . 
We point out that the behaviour of these amplitude and 
phase distributions is the same as the one obtained with the 
best solutions of the moment method although the values 
of the continuous parameters are slightly different. Again, 
we can exclude the (2, —2) solution while we cannot dis- 
criminate between a (2, —1) or (1, —1) mode. Since the pho- 
tometric mode identification in Paper I favours l\ = 2, we 
finally conclude that the wavenumbers of the main mode 



5 STELLAR PARAMETERS 

By means of the CORALIE echelle spectra we derive spec- 
troscopic estimates of the stellar parameters T e g and logp 
of the primary and tertiary by considering the presence of 
a B5 companion through the estimated flux ratio derived 
in Paper I. The used spectra were corrected for the orbital 
velocity obtained in Section 3 and some 20 spectra with the 
highest S/N ratio among the 86 CORALIE spectra at our 
disposal were selected to compute a mean spectrum with an 
increased S /N ratio and for which the pulsational variability 
of the primary is cancelled. We used five hydrogen lines (H a , 
Hp, H 7 , Hs, H.), six Hel lines (4026, 4387, 4471, 4713, 4992 
and 6678 A), two Sill doublets (4128-4131 and 5041-5056 
A), two Si III triplets (4553-4568-4575 and 4813-4819-4829 
A) and two additional Si III lines (4716 and 5740 A). 

These observed spectral lines were compared with theo- 
retical ones for both the primary and tertiary computed for 
spherically symmetric non-local thermodynamic equilibrium 
(NLTE) atmosphere models by me ans o f the FASTWIND 
code in the version of Puis et al.| (l2005h . The main input 
parameters to obtain the theoretical profile are T c g, logg, 
the microturbulent velocity and the particle number ratios 
n(He)/n(H), n(Si)/n(Si©). The gravity and temperature are 
most constrained by respectively the wings of the hydrogen 
lines and the ionization balance between Sill and III. 

The best overall fit to the observed line profiles after 
merging the theoretical profiles for the primary and tertiary 
occurs for the following parameters. For the primary and the 
tertiary, one has respectively T e g = 24000 ± 1000 K, logg = 
4.1 ± 0.1 and T cB = 18500 ± 1000 K, log g = 4.0 ± 0.1 with 
solar He and Si abundances. In Fig.HJwe show as examples a 
selection of observational hydrogen, helium and silicon line 
profiles together with their best merged line fits. Our results 



A stero seismology of 9 Oph: spectroscopy 7 



Ha 



1 .2 
1 .0 



0.4 




6540 6550 6560 6570 6580 6590 



4300 



4320 



4340 



4 560 



4 560 



Hel 4471 




Hel 4922 



4460 



4465 4470 4475 4480 




4915 4920 4925 4930 



.01 



1 .00 



0.99 r 



0.98 E_ 
4684 



Hell 4686 




Silll 4bb2-4b67-4574 




4685 



4686 



4687 



4688 



0.70 I 

4550 4555 4560 4565 4570 4575 4580 



Figure 6. Selection of observational hydrogen, helium and silicon line profiles together with their best fits obtained for T e g = 24000±1000 
K, log g = 4.1 ± 0.1 for the primary and T c jf = 18500 ± 1000 K, logg = 4.0 ± 0.1 for the tertiary and solar He and Si composition. 



of this spectroscopic determination of stellar parameters is 
in agreement with the photometric ones derived in Paper I. 



6 SUMMARY 

Recent thorough studies of a few (3 Cephei stars showed 
us that probing the internal structure of this kind of stars 
has become feasible. To this end, very intensive and long- 
term monitoring is required, as well as accurate determina- 
tion of many pulsation frequencies together with successful 
mode identification. For this reason, photometric and spec- 
troscopic multisite campaigns and long-term monitoring of 
several (3 Cephei pulsators are currently ongoing. 9 Ophiuchi 
is one of the studied targets for which the observational re- 
sults are presented in Paper I and in this paper. The ap- 
plication of asteroseismic techniques will be described in a 
following paper. 

We can summarize our findings as follows. The 1303 
colour photometric three-site data obtained during 77 nights 
of observation allowed the authors of Paper I to detect seven 
pulsation modes within a narrow frequency interval between 
7.116 and 7.973 c d _1 . Our less numerous data of 121 spectra 
were obtained during three years and allowed us to discover 
a low-mass spectroscopic companion of the B2 primary and 
to determine an orbital period of some 57 days. The pres- 
ence of a Speckle B5 companion was already known so that 



we can conclude that 8 Oph is a triple system. Moreover, we 
could recover three of the seven pulsation frequencies of the 
primary observed in photometry by means of a frequency 
analysis on the radial velocity but also from a two dimen- 
sional frequency search on the line profiles themselves. 

The photometric and spectroscopic mode identifications 
turned out to be complementary and fully consistent with 
each other. The colour photometry attributed unambigu- 
ously the frequency 7.4677 cd _1 to a radial mode, which is 
corroborated by spectroscopy. Two multiplets are also iden- 
tified in photometry: one complete 1=1 triplet, which is 
not observed in the spectroscopic data at our disposal, as 
well as three components of an t = 2 quintuplet, of which 
we detected two components in our spectroscopic data. Be- 
cause this quintuplet is not complete, there are still two 
possibilities for the m-values of the observed I = 2 modes. 
This ambiguity was lifted thanks to our spectroscopic mode 
identification, taking into account the photometric identifi- 
cation of the radial mode. By comparing the observed first 
three moments of the Si III 4553 A line as well as the ob- 
served amplitude and phase variations across this silicon line 
to the same theoretically computed quantities for several 
wavenumbers (I, m) and for a large grid of continuous pa- 
rameters, we identified the main mode with frequency 7.1160 
cd _1 as a non-axisymmetric (2, —1) mode and consequently 
we fixed the mode identifications of the whole quintuplet. 
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The position of the primary in the Hertzsprung-Russel 
diagram was also photometrically and spectroscopically de- 
termined, leading to compatible stellar parameter values. 

These observational results form a good starting point 
for in-depth seismic modelling of the star. In particular, since 
its frequency spectrum is so similar to the one of V836 Cen, 
we have the opportunity to test if the occurence of core over- 
shoo ting a nd non-rigid rotati on found for that star (Aerts 
et al. l200l Dupret et al. l2004l also applies to 6 Oph. 
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